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Abstract 

We present a method to compute non-perturbatively the renormalization constant 
of the scalar density for Ginsparg- Wilson fermions. It relies on chiral symmetry and is 
based on a matching of renormalization group invariant masses at fixed pseudoscalar 
meson mass, making use of results previously obtained by the ALPHA Collaboration for 
0(a)-improved Wilson fermions. Our approach is quite general and enables the renor- 
malization of scalar and pseudoscalar densities in lattice regularizations that preserve 
chiral symmetry and of fermion masses in any regularization. As an application we 
compute the non-perturbative factor which relates the renormalization group invariant 
quark condensate to its bare counterpart, obtained with overlap fermions at (5 = 5.85 
in the quenched approximation. 
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1 Introduction 



The understanding of the low-energy sector of QCD is one of the main goals of lattice 
simulations of the theory. In particular, the determination of the scalar quark-antiquark 
condensate associated with the spontaneous breaking of chiral symmetry has been the 
focus of many recent studies There are now good prospects for a reliable and pre- 

cise calculation of this quantity, owing to recent progress in two key areas. The first is 
the formulation of chiral symmetry on the lattice: it has recently become clear how the 
familiar consequences of chiral symmetry in the continuum can be preserved at non-zero 
lattice spacing, by employing lattice Dirac operators that satisfy the Ginsparg- Wilson 
relation || (for reviews of the subject see [pHIlf). In the context of spontaneous sym- 
metry breaking this makes it possible to extract the quark condensate in a conceptually 
clean manner through a suitable finite-size scaling analysis p, 12 1. Lattice results for 
the bare, subtracted condensate in infinite volume, — £ su b, have already been obtained 
in this way in the quenched approximation using Neuberger's (or overlap) opera- 

tor 12] as a realization of the Ginsparg- Wilson relation. In order to present an estimate 
for the condensate which can be used in phenomenological applications, the bare lattice 
result must, of course, be renormalized. 

The renormalization of matrix elements of composite operators defined on the lat- 
tice is the second key area where significant progress has been achieved (for a recent 
review see fli}] ). A theoretical framework to address the problem of (in general scale- 
dependent) renormalization of lattice operators in a completely non-perturbative man- 
ner has been developed. The main idea is the introduction of an intermediate renormal- 
ization scheme, such as the Regularization Independent (RI) (l]| and the Schrodinger 
functional (SF) [|l(| schemes. These formalisms have already been used successfully 
to address the renormalization of quark masses for Wilson [17-|l9|], staggered and 
Domain Wall [21| fermions. 

The main subject of this paper is the description of a method to compute non- 
perturbatively the multiplicative renormalization constant of the scalar density in fer- 
mionic regularizations based on the Ginsparg- Wilson relation. Our strategy relies on 
the matching of renormalization group invariant quark masses at fixed pseudoscalar me- 
son mass, and avoids the direct formulation of intermediate renormalization schemes, 
such as the SF, for Ginsparg- Wilson fermions. Instead, we make extensive use of the 
non-perturbative renormalization factor for quark masses, computed by the ALPHA 
Collaboration for 0(a)-improved Wilson fermions [jnj. The method is applicable to 
the renormalization of scalar and pseudoscalar densities for discretizations of the Dirac 
operator that preserve chiral symmetry at non-zero lattice spacing, and to the renor- 
malization of quark masses in any fermion discretization. As an example, we determine 
non-perturbatively the renormalization factor which links the bare scalar condensate 
obtained using the overlap operator to the renormalization group invariant condensate. 

The important issue of quenching will not be adressed here. While the methods de- 
scribed below can be applied without change to the full theory, the calculations discussed 
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are sufficiently costly with present computer and algorithm technology that unquenched 
calculations are not yet considered. Our main purpose is to demonstrate the feasibility 
of our strategy, and we restrict ourselves to a single value of the lattice spacing at which 
the bare subtracted condensate has been computed previously, corresponding to a bare 
coupling of (3 = 5.85. 

As we will describe below in some detail, our approach to the non-perturbative 
determination of the renormalization factor for the scalar density requires the calcula- 
tion of two-point functions in the pseudoscalar channel using overlap fermions. These 
correlation functions are themselves a rich source of information on low-energy QCD, 
and allow us to compute the condensate independently from the finite-size scaling anal- 
ysis of []|]. Furthermore, they serve to calculate the light quark masses as well as the 
pseudoscalar decay constant in the chiral limit. The detailed discussion of these results 



is deferred to a companion paper [22 



The outline of the remainder of this paper is as follows. In Section |2| we present 
our strategy for renormalizing the condensate. A discussion of cutoff effects associated 
with the use of an intermediate 0(a)-improved Wilson regularization and the intrinsic 
precision of our method is presented in Section [3|. Section || contains a description of our 
simulation results, using the overlap operator as well as 0(a)-improved Wilson fermions 
at (3 = 5.85. In Section |5| we present our results for the renormalization factors and 
compare them to one-loop perturbation theory. Finally, Section || contains a summary, 
including an estimate for the renormalized condensate at (3 = 5.85. Some details con- 
cerning improvement coefficients for 0(a)-improved Wilson fermions are described in 
Appendix |A[ 

2 Strategy 

In this section we describe the renormalization of the scalar condensate for overlap 
fermions. Our strategy relies on the fact that the renormalization constants for the 
scalar and pseudoscalar densities are identical and are equal to the inverse of the renor- 
malization factor for fermion masses, as can be shown using the chiral Ward Identities 
for Ginsparg- Wilson fermions (see, e.g. ref. 



Z P = Z S = (2.1) 

The non-perturbative renormalization of quark masses has been studied extensively on 
the lattice. In particular, the relation between the renormalization group invariant 
(RGI) mass and its counterpart in the intermediate SF scheme is known with an accu- 



racy of better than 2% in the continuum limit [24]. Furthermore, the non-perturbative 
matching between the SF scheme and 0(o)-improved Wilson fermions has been per- 
formed for a range of bare couplings |24fl. One possible strategy is then to repeat this 



second step for overlap fermions, by evaluating the normalization condition of ref. [24]. 



However, the direct implementation of the SF scheme for the overlap operator is not 
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straightforward, due to the inhomogeneous boundary conditions in the time direction, 
which are incompatible with the Ginsparg- Wilson relation. We have thus devised an 
alternative strategy, which allows us to exploit the previously obtained non-perturbative 
relations between quantities defined in the 0(a)-improved Wilson theory and their RGI 
counterparts. 

We begin by considering the renormalization constant Zm($o) which relates a bare 
quark mass, m(go), to the RGI mass M through 

M = Z M (g )m(g ), (2.2) 

where we have explicitly indicated the dependence on the bare coupling go . Note that 
we have not specified the fermionic discretization at this point. In the case of 0(a)- 
improved Wilson fermions we have 

M = Z^(g )m w (g ), (2.3) 

where m w is the current quark mass, and the superscript "w" reminds us that this 
defines Zm for this particular regularization. One can now write the ratio M/m{go) as 

M M m w (g' ) 



m(g ) m w (g' Q ) m(g ) 

7W [f J\ ( r m w)(ffp) , , 

— ^MOo ■ —, — 77 7"- 

(r m)(go) 

Here, g' Q is a value of the bare coupling which may differ from go, and in the last line 
we have also introduced the hadronic radius ro [^] to set the scale. The factor Z M has 
already been computed in the quenched approximation for a large range of couplings |24| . 
It is then clear that the relation between the RGI mass M and the bare mass m(go) 
in any regularization can simply be obtained by determining the values of (ro m) and 
(ro m w ) which reproduce a reference value rr re f of a chosen observable. A convenient 
choice is the pseudoscalar meson mass in units of tq, such that (fo^p) — 3?ref • 

Thus, 



■ W ( [Z^(g' )x(r m w )(gi )]. ' 



m( ff o) 1 v—w/^j {rom){go] 



(2.6) 

(r mp) 2 =x [c( . 



It is now important to realize that the combination Z^(g' ) x (ro r/7. w )(<7 ) is a renor- 
malized, dimensionless quantity. We can therefore define the universal factor Um in the 
continuum limit as 

U M = hm n {Z M (g' ) x (r m w )(g' )} (2.7) 



This completes our definition of the renormalization factor Zyi for any given fermionic 
discretization. By combining eqs. fl2.2j ), (^) and ( ^.7[) we obtain 

Z M (go) = U M ■ . , 2 (2.8) 

(rom) (r m P ) 2 =x rcf . 
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At this point all reference to the bare coupling g' Q and the use of 0(a)-improved Wilson 
fermions has disappeared, and the only part that retains an explicit dependence on 
the lattice regularization is the bare quark mass in units of Tq, (rom), at the reference 
point x re {. The only discretization errors that remain are those associated with the 
regularization for which Zm(9o) is considered. Estimates for Um in the continuum limit 
are easily obtained from published results employing 0(a)-improved Wilson fermions, 
which greatly facilitates the evaluation of Zm in any given regularization. We will return 
to this point in Section || 

Let us now consider a fermionic regularization which preserves chiral symmetry. 
In this case the renormalization factor Zm of eq. ( |2.8D serves not only to renormalize 
the quark mass, but also the scalar condensate. For concreteness we choose overlap 
fermions and assume that eq. ( |2.8| ) has been evaluated for m = m ov , where m ov denotes 
the bare mass in the massive overlap Dirac operator given below in eq. (4.1). [j The RGI 



condensate X is then obtained from the bare subtracted condensate computed using the 
overlap operator, X S ub(<?o)! through 

' S sub ( 50 ). (2.9) 



Zm(9o) 



The RGI condensate £ is a fully non-perturbative quantity. However, it is traditional 
to quote the value of the condensate in a perturbative scheme such as MS, at some 
reference scale \i. The matching of the RGI condensate to that defined in the MS 
scheme must necessarily be perturbative. The MS condensate Sj^g(^) is thus given via 



1 



%s(M) = v ( „_ ~\ S sub(5o), (2.10) 



Z m ( go ^) = ^d^Z M (g ). (2.11) 



where 



The numerical value for factor m^g(/i)/M was obtained in jllj] through numerical in- 
tegration of the perturbative renormalization group (RG) functions in the MS scheme. 
For instance, at the commonly used reference scale [i = 2 GeV the integration of the 
4-loop RG functions yields 



M 



0.72076, /x = 2GeV, (2.12) 



where the error due to the uncertainty in the quenched value of A^jg is 1.5%. The 
conversion to other reference scales is easily performed using the tabulated values of 
fn m (n)/M in Table 3 of ref. @]. 

3 Note that for Ginsparg- Wilson fermions the bare mass which appears in the Lagrangian is identical 
to the quark mass defined through the PC AC relation with the conserved axial current. 
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The relations between the bare and renormalized condensates of eqs. ( [2.9| ) and (|2.10[) 
are more conveniently expressed in terms of renormalization factors Z$ and Z$, which 
are related to Zm and Z m by 

Zs(9o) = ^r^, Z s (g ^) = —-^ -. (2.13) 

Our goal is to renormalize the bare subtracted condensate calculated in by com- 
puting Zs at (3 = Q/g^ = 5.85. Thus, in the remainder of this paper we will focus on 
the determination of (ro m ov ) for (5 = 5.85 at a suitably chosen reference point x re f. 
However, before describing the details of this calculation, we discuss some of the issues 
surrounding the determination of the universal factor Um- 

3 The universal factor Um , cutoff effects and overall precision 

The factor Um of eq. (2.7) is in fact the RGI quark mass (in units of ro) in the continuum 



limit for a degenerate, pseudoscalar meson with (Vqttzp) — %ref- 

In ref. @ RGI quark 

masses were computed for (rgmp) 2 = 1.5736 and 3.0. The first value corresponds to a 
degenerate pseudoscalar meson, which has the same mass as the kaon.0 This implies 
that Um = (M s +M)ro/2 for x re f = 1.5736, where M s is the RGI strange quark mass and 
M = \{M U + Md). On the other hand, choosing (romp) 2 = x re f = 3.0 corresponds to a 
degenerate meson with a quark mass roughly equal to that of the strange quark, such 
that Um ~ M s ro- In order to explore the systematics of our procedure more thoroughly 
we have considered a third reference point, x re f = 5.0, and the reasons for this choice 
are explained in Section |||. 



The results of ref. [19] are easily converted into estimates for Um a t £ re f — 1.5736, 
3.0 and 5.0: an estimate of (M s + M)ro in the continuum limit is given in eq. (5.2) 
of that paper. Furthermore, by evaluating eq. (5.1) for (romp) 2 = 3.0 in conjunction 
with the results in the last line of Table 2, one can infer the value of Um a t £ref = 3.0. 
Finally, using the results in Table 1, the procedure of [|HJ] can be repeated in order to 
determine Um for x re f = 5.0. Thus we obtain 

f 0.181(6), x re{ = 1.5736, 
U M = \ 0.349(9), x ref = 3.0, (3.1) 
[ 0.580(12), rc ref = 5.0. 

Therefore, what is left to do in order to obtain a fully non-perturbative mass renorma- 
lization constant Zm in any given regularization is to calculate (ro m) in that scheme 
for one (or all) of these reference values. 

Before undertaking this calculation, we wish to discuss the issue of discretization 
errors in our renormalization condition, eq. (|2.8|) . To this end we note that a valid 
definition of Zm is also provided by the right-hand side of eq. ([[!]) evaluated at non- 
zero g' . The two expressions differ by 0(a 2 ) discretization errors associated with the 
4 Using m K = 495 MeV and r Q = 0.5 fm gives (r m K ) 2 = 1.5736. 



5 



0.6 



0.5 



£ 0.4 



o 



s 

IS1 



0.3 



0.2 



0.1 




ref 

_1 I L 



1.5736 

I I I L 



_1 I L_ 



£11 

J L 







0.02 0.04 
(a/r ) 2 



0.06 



Figure l: Cutoff effects in Z M (ro m w ) for x re f = 1.5736, 3.0 and 5.0. The continuum 
estimates for Um of eq. ( \3. \ ) are represented by the error band (dotted lines). The 
corresponding linear fits in (a/ro) 2 to data computed for f3 > 6.1 are shown as the solid 
lines. Extrapolations using data at all five values of the lattice spacings are represented 
by the dashed lines. 



intermediate Wilson regularization. In eq. (|2.6| ) these errors are present, while they are 
not in eq. fl2.8|) , owing to the continuum extrapolation in the definition of Um- One 
might think that such 0(a 2 ) artefacts - which are formally of the same or even higher 
order than those of most other lattice regularizations - are unimportant. We are now 
going to show that this is not the case. To this end we have evaluated 



z m(9o) x (rom w )(g' 0/ 



(romp) 2 =i ro f, 



x re f = 1.5736, 3.0, 5.0 



(3.2) 



for bare couplings g' corresponding to (3 = 5.85,6.0,6.1,6.2 and 6.45. Data for (3 > 6.0 
were taken from ref. [19|. The results for (r m w ) at /3 = 5.85 were obtained as described 
in Section || below. The value for Z M at /? = 5.85 was estimated by extrapolating the 
parameterization of eq. (6.10) in [24 and doubling the error. 

The approach to the continuum limit is shown in Fig. [l|. The first observation 
is that lattice artefacts in Z M (rom w ) appear to be consistent with the expected a 2 
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behaviour - at least for x re f = 1.5736 and 3.0: performing linear fits in (a/ro) 2 to data 
points at all five values of the lattice spacing produce good x 2 /dof. Furthermore, as 
can be seen from the figure, the results in the continuum limit are compatible with the 
estimates for Um in eq. (|3.1| ). The latter were obtained by excluding all points below 
f3 = 6.1 from the continuum extrapolations, as a safeguard against potentially large 
cutoff effects of higher order. As can be seen from the figure, such discretization effects 
appear to be larger at x re f = 5.0, which certainly justifies the exclusion of those data 
points obtained on the two coarser lattices. 

The deviation of Z M (rQm w ) from its continuum value amounts to about 20% at 
j3 = 6.0 and 40% at j3 = 5.85. Thus, if Zm were evaluated using eq. (^) with g' = go 
instead of eq. ( |2.8D , the result at (3 = 5.85 would be 40% smaller, due entirely to cutoff 
effects of order a 2 in the intermediate Wilson regularization. 

This discussion underlines that it is important to use the continuum result Um 
in the renormalization condition. It guarantees that the approach to the continuum 
limit of quantities for which these renormalization constants are used is not obscured 
by lattice artefacts of 0(a 2 ), introduced by using 0(o)-improved Wilson fermions as an 
intermediate regularization. 

As can be inferred from eq. (|3.1|), the factor Um is known with a precision of about 
3%. At present this presents a lower bound on the accuracy in the determination of 
renormalization factors according to our proposal. One might argue that more precise 
estimates for Um could be obtained by including the data with (3 < 6.0 in the continuum 
extrapolation. We have still decided against this, in order to be certain about excluding 
effects from higher orders in the lattice spacing. Also, given the high cost of simulations 
employing discretizations such as overlap fermions, it will be some time before the lower 
bound of 3% will be regarded as a real limitation. If one wants to improve the accuracy 
in the determination of Um, it would be preferable to add data points at smaller lattice 
spacings, which should be possible with a relatively modest amount of computer time. 

4 Lattice calculation of pseudoscalar masses 

We now describe the details of our calculation of pseudoscalar two-point functions using 
the overlap operator as well as 0(a)-improved Wilson fermions at (3 = 5.85. Although 
our strategy only requires published results for Wilson fermions, the additional Wilson 
calculations performed for different lattice sizes provide us with valuable information 
about finite-size effects at a much lower cost than with overlap fermions. 

4.1 Simulation details 

For the massive overlap operator we use the following definition 

^=(l-^^y)^ 0) +-o V , (4.1) 



7 



aD { $ = (1 + s) ( 1 - A/VaTa) , A=l + s-aD w , (4.2) 



where 

and Dw is the standard Wilson-Dirac operator. The naive continuum limit of Dn is the 
canonically normalized Dirac operator with a bare quark mass m ov . The parameter s 
was fixed at s = 0.6, which is close to the value where the localization of the operator 
was found to be optimal at this j3 value pq|. The numerical implementation of the 



inverse square root in eq. (44) was performed using a Chebyshev approximation, and 
for further details we refer to our earlier work ||. Here we only mention that we have 
employed a multi-mass solver [27], which allows for a simultaneous calculation of the 



quark propagators for seven values of the bare mass ranging from am ov = 0.047 — 0.188. 

In order to compute pseudoscalar meson masses and current quark masses for 0(a)- 
improved Wilson fermions we have followed the same procedure as in ref. |2£|. In 
particular, we have employed Schrodinger functional boundary conditions on lattice 
sizes L/a = 8, 10, 12 and 16, with T/a = 24 fixed, in order to study finite-size effects 
For the improvement coefficients c sw and ca, which appear in the definition of the 
improved action and axial current, respectively, we have chosen 

c sw = 1.909, c A = -0.144. (4.3) 



The value for c sw was obtained from the interpolating formula eq. (4) of ref. pS|] . To our 
knowledge a non-perturbative result for the coefficient ca has so far not been published 
for f3 < 6.0. The choice in eq. (|4.3|) is based on a non-perturbative calculation using the 
SF, which is described in detail in Appendix [A]. 

While it can be shown that there are no exceptional configurations for the overlap 
operator at non-zero quark mass [13,|30, 31 1, this is not the case for Wilson fermions. 



Indeed, when working below (3 = 6.0 the incidence of exceptional configurations may be 
so high - especially for quark masses below that of the strange quark - so as to make a 
reliable calculation of quark and meson masses impossible. 

For our calculations of quark propagators using Wilson fermions we have cho- 
sen relatively heavy quarks, with unrenormalized current quark masses in the range 
80 — 120 MeV. Only on the smaller volumes of L/a = 8, 10 did we push to lighter 
quarks, corresponding to 65 MeV. We checked against the occurrence of exceptional 
configurations by plotting the Monte Carlo history of the correlation functions of the 
axial current and pseudoscalar density evaluated at xq = T/2. Exceptional configura- 
tions manifest themselves as isolated peaks (or dips) whose heights exceed the typical 
statistical fluctuations by several orders of magnitude. Although one candidate each was 
detected for L/a = 8 and 10, the observed fluctuations were not deemed large enough 
to justify the exclusion of these configurations from the statistical ensembles. Hence, 
for the case of Wilson fermions, all hadron and quark masses have been evaluated using 
the full statistics on all lattices. 



3 As an orientation for the reader we add that L/a = 12 corresponds to L ~ 1.5 fm. 
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L/a 


configs. 


K 


am w 


amp 


(amp) 2 


8 


640 


0.13150 
0.13200 
0.13250 
0.13300 


0.07309(37) 
0.06106(40) 
0.04893(47) 
0.03685(52) 


0.5780(45) 
0.5265(49) 
0.4707(54) 
0.4085(60) 


0.3342(52) 
0.2772(51) 
0.2216(50) 
0.1668(49) 


10 


512 


0.13150 
0.13200 
0.13250 
0.13300 


0.07362(27) 
0.06174(28) 
0.04991(30) 
0.03805(34) 


0.6002(29) 
0.5496(31) 
0.4948(35) 
0.4338(40) 


0.3602(35) 
0.3021(35) 
0.2449(35) 
0.1882(35) 


12 


256 


0.13150 
0.13200 
0.13250 


0.07436(30) 
0.06249(32) 
0.05067(35) 


0.6005(30) 
0.5503(33) 
0.4961(36) 


0.3606(36) 
0.3028(36) 
0.2461(36) 


16 


150 


0.13150 
0.13200 
0.13250 


0.07458(24) 
0.06277(25) 
0.05102(27) 


0.6065(21) 
0.5573(23) 
0.5045(26) 


0.3678(26) 
0.3106(26) 
0.2545(26) 



Table l: Lattice sizes, statistics, pseudoscalar and current quark masses computed using 
0(a)-improved Wilson fermions. 

4.2 Results for pseudoscalar masses 

The results for the pseudoscalar and bare current quark masses computed for 0(a)- 
improved Wilson fermions are presented in Table [l]. The masses were extracted follow- 
ing the procedure described in detail in ref. |28[| . In particular, the estimates for the 
pseudoscalar masses listed in the table were obtained by averaging the effective masses 
computed from the correlation function of the improved axial current. In accordance 
with Table 1 of |f2£f| , and using r$/a = 4.067 ± 0.014 [^] we have chosen the time 
window 11 < xq/o < 15 for the averaging procedure. Estimates for the bare current 
quark mass were obtained in a similar manner, by averaging the results over the interval 
10 < xq/o < 16. This choice of time window coincides with a clear plateau observed for 
this quantity. 

By comparing the results in Table [l] obtained on different lattice sizes, one observes 
that the time interval 11 < xo/a < 15 suggested in [^] produces consistent results for 
pseudoscalar masses if L/a > 10. By contrast, using a spatial lattice size of L/a = 8 
(corresponding to L ~ 1 fm) leads to finite-size effects at the level of 4% (4 standard 
deviations). Choosing a shorter time window on L/a = 8, such as 11 < xo/a < 13, 
brings the results closer to those on the larger volumes, but a 2% effect (1.5 standard 
deviations) remains. 

On the other hand, for L/a > 10 (i.e. L > 1.2 fm) finite-size effects in the pseu- 
doscalar mass are of the order of 1.5% or less. In order to exclude large finite-size 
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am ov amp (amp) 2 



0.047 
0.063 
0.078 
0.097 
0.125 
0.161 
0.188 



0.311(31) 
0.357(22) 
0.398(17) 
0.444(14) 
0.506(12) 
0.579(11) 
0.630(10) 



0.096(20) 
0.128(16) 
0.159(14) 
0.197(13) 
0.256(13) 
0.335(12) 
0.397(12) 



Table 2: Pseudoscalar masses extracted from single-mass fits in the interval 6 < x^/a < 
12 at several values of the bare mass m ov . 

effects for pseudoscalar masses computed using the overlap operator, whilst keeping the 
computational overheads small, we have thus decided to work with L/a = 10. However, 
the above analysis of finite-size effects only applies to pseudoscalar meson masses with 
amp 0.5, a, range which, lies above the reference points x re f — 1.5736 and 3.0. This 
was our main reason to consider the reference point x re f = 5.0 in addition: at (3 = 5.85 
this choice corresponds to amp ~ 0.55, a value for which the absence of significant 
finite- volume effects has been confirmed. 

In Table ^, we present the results for pseudoscalar masses computed using the 
overlap operator on an ensemble of 50 configurations. The chosen lattice size was 
L/a = 10, T/a = 24, with periodic boundary conditions in all space-time directions. 
After averaging the correlation functions over the forward and backward directions in 
euclidean time, we extracted the masses from single-cosh fits to the correlation functions 
in the interval 6 < xq/o < 12. The quoted statistical errors were obtained from a 
jackknife procedure. We have classified configurations according to their topological 
index, distinguishing between topological (having a non-zero index) and non-top ological 
configurations and determined the average pseudoscalar masses restricted to either class. 
Only very near the chiral limit does one expect these quantities to differ. In the range 
of quark masses considered we did not detect any difference between the two classes, so 
that we could safely include all topologies in the average. 

In Fig. ^ we show the results for the behaviour of (amp) 2 as a function of am ov 
and am w . In both cases the results are perfectly compatible with the linear behaviour 
expected in lowest order Chiral Perturbation Theory. As an illustration we have plotted 
the results from a chiral fit to our data, using the linear parameterization 



thus allowing for a non-zero intercept of (amp) 2 at vanishing quark mass. The results 
of these fits are: 




(4.4) 



A ov = -0.009(18) 
A w = 0.0039(39) 



B ov = 2.139(80) 
B w = 4.836(44), 



(4.5) 
(4.6) 
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Figure 2: irip versus m ov (open circles) and m w (filled circles) respectively for overlap 
and O (a) -improved Wilson fermions. The solid lines are the fits of eqs. jj.^j) and 



where we have used all seven data points in the overlap case and all four quark masses 
in the Wilson set. Note that in both cases the intercepts are perfectly compatible with 
zero within errors. 

Pseudoscalar propagators for overlap fermions have been studied before in for 
j3 = 5.7,5.85 and 6.0, thus enabling a direct comparison with our results.^]. Since the 
authors of [33] do not present numbers for (amp) 2 one has to infer its quark mass 
behaviour by reading the data off their plots. From this fairly crude comparison we 
conclude that the slope parameter B ov is in rough agreement with our findings. Unlike 
the authors of [33] we did not attempt to model the quark mass behaviour including 



quenched chiral logarithms. Judging from the quality of our fits to eq. (4^) we do not 
expect a significant deviation from leading order Chiral Perturbation Theory at our level 
of statistics in the range of quark masses we consider. In order to test for the presence 
of quenched chiral logarithms, it may be necessary to go to much smaller masses than 
our lightest quark, whose mass is roughly half as large as that of the strange quark. 



5 Results and discussion 



6 Note that the bare mass in |33| differs from m 
a factor (1 + s). 



To leading order in a the two definitions differ by 
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5.1 Non-perturbative result for Zs 

We have now the necessary results to compute the renormalization factor Zs for the 
RGI condensate using eqs. ( |2,8| ) and ( 2 . 1 3| ) . Our task is to determine (rom ov ) at the 



reference value x re f , which could in principle be achieved using the results of the chiral 
fit, eq. d4.5|) . However, we prefer to perform local interpolations of the quark mass 
m ov to the reference point, which avoids any assumption about the mass behaviour of 
(amp) 2 at or very near the chiral limit. 

To this end we have interpolated (rom ov ) to our chosen values of x re f by using the 
three nearest data points for (romp) 2 . We obtain 

f 0.190(44), x ref = 1.5736 
(r m ov h ^ =\ 0.363(26), x ref = 3.0 (5.1) 
(r ° mp) =Xrcf { 0.594(23), x ref = 5.0. 

We add that these results are entirely consistent with interpolations using only the two 
neighbouring points. Given the almost perfect linearity of the quark mass dependence 
of (romp) 2 , it is not surprising that interpolations using all seven data points are also 
consistent. Hence we regard the uncertainties in our results of eq. (|5.1| ) as conservative 
estimates. 

By combining the results for (tq m ov ) with the factor Um we obtain 

f 1.05(25), x re{ = 1.5736 
Z s = J 1.04(8), x re{ = 3.0 = 5.85. (5.2) 

[ 1.02(4), x ref = 5.0, 

The results at the three values of 

%xef £11*6 consistent within even the smallest of the 
statistical errors. This indicates that our renormalization condition may be applied 
over a fairly large range of quark masses without introducing significant discretization 
errors due to working at non-zero quark mass. It also indicates that finite- volume effects, 
even at the lightest reference point, appear to be small. Owing to the modest statistics 
in our simulations, the errors in Zo, in eq. ( |5.2| ) are clearly dominated by the limited 
accuracy of the estimates for (ro m ov ) at the reference points, which amounts to 24%, 
8% and 4% at x re f = 1.5736, 3.0 and 5.0, respectively. This is larger than the relative 
error in Um of about 3%. However, with better statistics, it should be easy to obtain 
more accurate estimates for (rom ov ) and thus Zs- 

The dependence on the value of x re f is quite weak and well covered by the statistical 
uncertainty. We quote the result for x re f = 3.0 as our best estimate, since this reference 
point is a good compromise between being chiral enough so as not to introduce significant 
discretization errors, and massive enough to guarantee negligible finite-size effects. Thus 
we obtain 

Z S =1.04 ±0.08 
Z s (2GeV) =1.44 ±0.11, P ~ 5 ' 85 ' (5 ' 3) 

These numbers are the main result of this paper. All systematics associated with the in- 
termediate Wilson regularization have been eliminated by the continuum extrapolation 
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in the definition of the universal factor Um ■ As discussed in Section |], we have verified 
that finite- volume effects are under control. Discretization errors of order (am) 2 in our 
determination of Z$ seem to be small, as indicated by the consistency of our results at 
the three reference values. 



5.2 Comparison with perturbation theory 



In ref . [ 23 1 the perturbative renormalization of quark bilinears was studied for the overlap 
operator. The one-loop expression for Zg reads 



(90, iA 



l + 5o 



where the one-loop coefficient depends on the parameter s in the definition of the 
overlap operator. For our value of s = 0.6 one finds 



7^2 ln ( a ^) + z s ] 



+ O(g A ), 



(5.4) 



0.107074. 



(5.5) 



It is well known that perturbation theory in the bare coupling go is not very convergent. 
As a consequence it has become customary to consider "mean-field improved" estimates 
for perturbative renormalization factors [34, 35]. Another proposal is based on the 
resummation of "cactus" diagrams [p3|| . Here, in the spirit of [35], we propose the 
following mean-field improved expression: 



l + s 
1 + 3 



1 + g 2 



^2 + 4 1} + 4 1} v , + s 



(5.6) 



where 3 = 3 + (s — 3)/uo- We take Uq to be the average plaquette in infinite volume, and 
Ug 1 ^ = —1/12 is the one- loop coefficient in its perturbative expansion. At this order, the 
choice of coupling, g, is ambiguous. For consistency, we work with the coupling used 
in obtaining the renormalization group factor in eq. ( 2.12| ), i.e. g 2 = 2.5432. We can 
now evaluate eq. ( |5.6| ) using lattice data for the average plaquette and by setting the 
lattice spacing at /3 = 5.85 with the help of the interpolating formula for the hadronic 



radius ro/a [32|. For u 



0.575 the results for Z™ f and Z™ f are 



Z^ f (2GeV) 

Z S 



1.26 
0.91, 



j3 = 5.85. 



(5.7) 



Comparing these mean-field improved perturbative estimates to the non-perturbative 
results of eq. (|5.3[), one finds that the former are about 12% smaller. This is a significant 



improvement on the results of bare perturbation theory, where Zg t (2GeV) = 1.12 as 



given by eq. (| 
of eq. (l5lf). 



, which are more than 20% smaller than the non-perturbative results 
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6 Summary and conclusions 



We have proposed and tested a method to non-perturbatively renormalize scalar and 
pseudoscalar densities for fermionic discretizations which preserve chiral symmetry at 
non-zero lattice spacing. This method also provides a means to renormalize quark 
masses non-perturbatively in any lattice regularization. Given the cost of simulations 
employing Ginsparg- Wilson fermions and the incompatibility of Schrodinger functional 
boundary conditions with the Ginsparg- Wilson equation, we found it advantageous to 
proceed through 0(a)-improved Wilson fermions, where the non-perturbative renorma- 
lization of currents and densities has already been studied extensively. The sought-after 
renormalization constants are then obtained through a matching of RGI quark masses at 
fixed pseudoscalar meson mass. As evident in eq. (|2.8|) , all reference to Wilson fermions 
drops out, owing to the universal factor Um defined and evaluated in the continuum 
limit. Furthermore, the overhead for implementing the matching condition should be 
negligible in most cases, since it relies on quantities that are commonly computed in 
phenomenological studies where the renormalization constants are likely to be used. The 
idea to use an intermediate lattice regularization which is relatively cheap to implement 
may prove useful for the computation of other renormalization constants in fermionic 
discretizations that are numerically much more demanding. 

As an application, we have computed the renormalization constants which are re- 
quired for the quark condensate obtained using overlap fermions in the quenched approx- 
imation. Our results at f3 = 5.85 are listed in eq. (|5.3| ) and can now be combined with the 
result for the subtracted bare condensate determined in ||, i.e. a 3 £ su b = 0.00323(37) []. 
In units of ro we obtain 



rjjS = 0.226(26)(2)(17) 
r 3 S M g(2GeV) = 0.313(36)(3)(23), 



(3 = 5.85, (6.1) 



where the first error is due to the statistical uncertainty in a 3 E su b, the second corre- 
sponds to the error in ro/a, and the third arises from the error in the renormalization 
factors Zs and Z$, respectively. Combining all but the error associated with ro/a in 
quadrature and using ro = 0.5 fm we find 



(240 ± HMeV) 3 x 



v 1605 MeV J 

E Hg (2GeV) =(268±12 MeV)S X fe^aY 



(3 = 5.85. (6.2) 



These results are still subject to discretization errors of 0(a 2 ). A detailed discussion 
of lattice artefacts and other systematic errors - including the scale ambiguity in the 

7 In ref. ]Bj the massive overlap operator was not O(a) improved. The above result for a 3 S su b was 
obtained for the O(a) improved definition in eq. ( ff.f| ), which amounts to a redefinition of the tree level 
mass. The difference from the result quoted in M is less than 1%. 
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Figure 3: The chiral extrapolation of ca using the three heaviest quark masses. Data for 
the standard and improved definitions of the lattice derivative are denoted by the open 
and filled symbols, respectively. 



quenched approximation - as well as a comparison of these results with those obtained 



through other approaches is deferred to our companion paper [22]. 
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A Determination of ca at f3 = 5.85 

While ct non-perturbcitive value for the improvement coefficient c sw is available at (3 = 
5.85 [^9), the coefficient ca of the improved axial current, which is required for the 
current quark mass, is not known for (5 < 6.0. 

As mentioned in Sect. numerical simulations using Wilson fermions are ham- 
pered by the occurrence of exceptional configurations associated with unphysical zero 
modes of the Wilson-Dirac operator. This problem is further exacerbated by working 
at a ^ 0.1 fm (i.e. (3 % 6.0) and small quark masses. However, if one is willing to relax 
the requirement that the improvement conditions for c sw and ca be evaluated at or very 



near the chiral limit [36|, the occurrence of exceptional configurations may be sufficiently 
suppressed. Thus, an extension of the determination of improvement coefficients to the 
regime where (3 < 6.0 may be possible. 
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Figure 4: The result for ca at (3 = 5.85 (filled circle) compared with the results from 
ref. Jffg/ (open circles). The solid curve is the interpolating formula valid for 6.0 < (3. 
The dashed curve denotes perturbation theory in the bare coupling, whereas the crosses 
represent the estimates in mean-field improved perturbation theory. 



In order to determine ca at (3 = 5.85 we have chosen c sw = 1.909 [29] and followed 



the strategy of ref. [36]. It has been observed pTj], though, that the original improvement 



condition for ca used in |36| suffers from a loss of numerical accuracy towards larger 
values of the bare coupling, corresponding to (3 ~ 6.0. A modified, but closely related 



improvement condition for ca has been proposed and tested 37]. It was shown that 
the alternative condition's numerical sensitivity does not deteriorate at large couplings. 
Furthermore, it gives consistent results compared with the original condition, thereby 
confirming the determination of ca in ref. |36| . 

We have evaluated both the original and modified improvement condition at (3 = 
5.85 on a 8 3 • 16 lattice. We used the same four values of the hopping parameter listed 
in Table |l|. With the procedure outlined in Sect. |] to check against the occurrence of 
exceptional configurations, we eliminated two configurations from the original ensemble 
of 2560. It was found that the alternative improvement condition gave a stable signal for 
ca at the three heaviest quark masses. By contrast, the original improvement condition 



used in [36] performed so badly at (3 = 5.85 that it could not be used to determine ca- 
Our estimates for ca obtained at non-zero quark masses shows some dependence 
on the quark mass, and therefore our final result is obtained through a chiral extrap- 
olation as shown in Fig. ||. We have also evaluated the improvement condition using 
a higher-order lattice derivative [Q, which was found to have a significant impact on 
the determination of some improvement coefficients The corresponding results and 
extrapolation are also shown in Fig. ||. Both definitions of the lattice derivative give 
entirely consistent results in the chiral limit, and as our final result we quote 

c A = -0.144 ±0.011, = 5.85. (A.l) 
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m 



This value is shown together with the previous determination of ca for /3 > 6.0 |3| 
Fig. | 

Finally we note that we have also estimated ca using the method proposed in 
11) 



ref. [40 (see also [41,[42|]). We have confirmed the observation of [}42(| , namely that the 
accuracy of the method is limited by the relatively small range in xq in which a stable 
signal is obtained. On our 8 3 • 16 lattice the method of |^0| produces results which are 
entirely consistent with the estimate in eq. ( A.l). 
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